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It is shown that the dynamics of domain walls in Left-Right symmetric models, separating respec- 
tive regions of unbroken SU{2)l and S'I7(2)_r in the early universe, can give rise to baryogenesis via 
leptogenesis. Neutrinos have a spatially varying complex mass matrix due to CP- violating scalar 
, condensates in the domain wall. The motion of the wall through the plasma generates a flux of 

(N . lepton number across the wall which is converted to a lepton asymmetry by helicity-flipping scat- 

terings. Subsequent processing of the lepton excess by sphalerons results in the observed baryon 
asymmetry, for a range of parameters in Left-Right symmetric models. 
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I. INTRODUCTION 



o 

\ Explaining the observed baryon asymmetry of the Universe within the framework of gauge theories and the standard 
' Big Bang cosmology remains an open problem. The study has resulted in a deeper understanding of nonperturbative 
rS I phenomena at finite temperature in gauge theories including supersymmetric theories. Many of the particle physics 
JL models and scenarios considered so far seem to require unnatural extensions or very special choices of parameters for 
^ successful baryogenesis; prime among these are the standard model (SM) and its minimal supersymmetric extension 
(~| (MSSM), using a first order phase transition ||^, Among the alternative proposals are those which rely on the 

L| ' presence of topological defects, viz., domain walls and cosmic strings ^. The latter are generic to many gauge 
theories. What makes them especially suited for baryogenesis is their nonthermal nature soon after their formation. 
Unlike the need for a first order phase transition which sets severe limitations on the couplings and particle content 
5-H ' of the model, existence of defects relies only on topological features of the vacuum manifolds and permits nonthermal 
effects without fine tuning. 

Many special features arise when studying cosmological consequences of topological defects in any given gauge 
model. The Left-Right symmetric (L-R) model was studied in the context of conventional baryogenesis mechanisms 
in 10, ||] and in the context of domain wall mediated mechanism in |^ . A detailed study of the possible defects existing 
in the L-R model was made in [p^ . R was argued that the domain wall configurations implied by the symmetry 
breaking pattern present possibilities for baryogenesis. In this paper we study the interaction of neutrinos which 
derive Majorana masses from the scalar condensate which constitutes the domain wall. Many of the broad features 
encountered, e.g., asymmetric reflection and transmission of fermions from moving domain walls, have appeared in 
the study of electroweak baryogenesis. In the diffusion-enhanced scenario driven by thick walls, the asymmetry 
diffusing in front of the wall is equlibrated by high temperature sphalerons. In our mechanism this is replaced by 
helicity flipping interactions in front of the wall which arise from the scalar condensate imparting a Majorana mass to 
the fermions. Our parametric answer for the unprocessed Lepton asymmetry produced in this mechanism is therefore 
dominated by the Majorana mass parameter / in eq. (^). The scalar condensate is absent behind the wall and 
therefore the asymmetry that has streamed through persists. 

At the completion of the L — R symmetry breaking transition, a particular hypercharge Y — — \{B — L) \s 
demonstrated to be spontaneously generated in the form of left-handed neutrinos. Due to the high-temperature 
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electroweak sphalerons, which set B + L = 0, this wiU be converted into an asymmetry of baryon minus lepton 
number {B — L). The baryon asymmetry thus generated arises in addition to that from the weU- known leptogenesis 
pj] , |l2j mechanism due to Majorana neutrinos. However, the two mechanisms constrain the Left-Right model rather 
differently. The usual mechanism requires the Zji mass to be larger than the heavy neutrino mass |l^, 
The present mechanism constrains the parameters of the Higgs sector for adequate CP violation, and the Majorana 
Yukawa couplings as already pointed out. Our main result is the identification of broad ranges of these parameters that 
ensure sufficient lepton asymmetry. Further, the subsequent evolution of this asymmetry must successfully produce 
the abserved baryon asymmetry. This requirement can be used to constrain the temperature scale of the L — R phase 
transition, eq. (^ ), or alternatively, the light neutrino mass, eq. (pO|). 

In [p^ and [|16| the possibility of B — L generated by any mechanism being neutralised due to presence of heavy 
majorana neutrino was considered. The bound first obtained in is 20il/„ > ^JTb-lMpi with M„ the mass scale 
of the heavy majorana neutrino and Tb-l the temperature at which B ^ L originates. This is derived from lepton 
number depletion due to heavy neutrino mediated scattering processes and assumes Tb^l > M„. It was argued in 
that the requirement that the heavy neutrino decays occur only in out of equilibrium conditions places a more 
stringent bound. Using the see-saw relation, it requires 

(^1/2 

< TO* = AngV^^ lO'^eV. (1) 

where Gj^ is the Newton constant and G^? is the Fermi constant. Current neutrino data easily suggest a larger 
neutrino mass. In this case it is argued that fl^ one needs 
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Tb-l < M„ - — — 10"GeV (2) 
2 \m,yj 

These considerations generically need a low scale for B — L creation. Detailed investigations |T7| of Leptogenesis 
scenrios, including lepton generation mixing, show that in several specific unified models this can be achieved in the 
context of conventional Leptogenesis. The present mechanism has the potential of meeting the requirement of low 
scale B — L generation in a natural way, although detailed investigations remain to be carried out. We shall return 



to this point in sec. |VI. 

The paper is organised as follows. Section || introduces the main features of the Left-Right symmetric model. 



synthesizing the conventions used by previous authors with the ones we follow. Section III discusses the microscopic 



mechanism of lepton number violation in scattering near the domain wall. It outlines the method that can be used 



for detailed study of lepton number creation in this context. Section IV demonstrates the existence of the conditions 
required for lepton number creation, in particular the CP-violating nature of the wall profiles. Section ^ presents 
a simplified version of the full theory to be studied and numerical results justifying the general conclusions of the 



previous section. Section VI discusses the implications to cosmology. Overall conclusions are presented in the last 
section. 



II. THE LEFT-RIGHT SYMMETRIC MODEL 



For the purpose of model building, Left-Right symmetry is a broad category, with several possible implementations. 
In this paper we shall adopt its more popular version which is desribed below. From the point of view of our mechanism, 
the discrete symmetry under exchange of the SU{2)b. field content of the model with SU{2)l field content is crucial. 
The breakdown of this symmetry gives rise to domain walls whose field configuration we study in detail. The most 
elegant version of the model consists of identical values of the two SU{2) gauge couplings in addition to a strict 
equality of certain scalar couplings in the Higgs potential. This may seem like an artificial requirement, considering 
that the two semisimple groups are independent, and there are no dynamical hints why they must be exactly same to 
begin with. More importantly, if this requirement is imposed an unpleasant feature arises in the context of cosmology. 
Breakdown of an exact discrete symmetry gives rise to stable domain walls and unless some mechanism removes them, 
they quickly come to dominate the energy density of the Universe, contrary to observations. Thus departure from 
exact symmetry is in any case a phenomenologically desirable feature. Happily, the mechanism being proposed here 



3 



[VA 



works well so long as the departure from exact symmetry is s mall s o that domain walls indeed form as transient 
constructs. A quantitiative discussion of this is taken up in sec. 

We now recapitulate the minimal SU{2)l SU{2)ji (g) U{1)b- 
scale vr, 



_i model ^8, 19 1. Parity is restored above an energy 
taken to be much higher than the electroweak scale, by introducing the SU{2)ii gauge symmetry which 
breaks at vr. Accordingly, a right-handed heavy neutrino species is added to each generation, and the gauge bosons 
consist of two triplets = (3, 1, 0), = (1, 3, 0) and a singlet = (1, 1,0). A Left-Right symmetric assignment 
of gauge SU(2) charges to the fermions shows that the new hypercharge needed to obtain the usual electric charge 
correctly is exactly B — L. It is appealing that in this model the weak hypercharge is related to known conserved 
charges. 

The electric charge formula now assumes a Left-Right symmetric form 



B-L 



(3) 



where and are the weak isospin represented by r'^/2, and is the Pauli matrix. 

The Higgs sector of the model is dictated by two considerations: the pattern of symmetry breaking and the small 
masses of the known neutrinos via the seesaw mechanism. The minimal set to achieve these goals is 



A, = 




(0,1,2) 



(4) 



where the electric charge assignment of the component fields has been displayed and the representation with respect 
to the gauge group is given in standard notation. 

The minimal form of the Higgs potential needed to fulfill the main phenomenological requirements can be found in 
pgf . This is however not the most general form. In ||2^ as well as the possibility of spontaneous CP violation was 
considered. In this case the couplings are chosen to be real, yet the translation invariant minimum of the potential 
occurs for complex VEV's (vacuum expectation values). The absence of explicit CP-violating couplings makes it 
easier to accommodate phenomenological constraints on CP violation. In the cosmological context in which we treat 
this theory, this motivation is not as compelling. Nevertheless, the same simplifying assumption will be made here. 

Consider the potential parameterized as [Eol 



with 



V = V^ + Va + V,j,a , 

14 = - filTrclxj)'' - Ai^(Tr^0t + Tr00^) 

+ Ai(Tr(/.0t)2 + A2{(Tr^(^t)2 ^ (Tr^t^)2| 

+ A3(Tr0(/)^)(Tr0^0) -f A4{Tr0V(Tr<?;.t0-|-Tr0l'0)} , 



(5) 



(6) 



Va= - A*iTr(At A, + At A,) + p, [(Tr(At A,)) + (Tr(At A.) 

TrAt AtTrA.A, + TrA^At TrA^^A J + pa [TrfAiA^) TrfA^ A, 



P2 



P4 TrA^Ai.TrAtA^-HTrAiAiTrAHA^ 



(7) 
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Ki,A = «iTr0(/)t(TrA,At +TrA,At) 

+ a2 {Tr(^t^)Tr(A,At ) + Tr(0(/.t)Tr(A,At)} 

+ al {Tr(00t)Tr(AKAt ) + Tr(^t0)Tr(A,At )} 
+ as {Tr(00tA,At) + Tr(0t0A,At )} 

+ /3i {Tr(0A,0tAt) + Tr(0tA,0At )} + {Tr(0A„0tAt) + Tr(^t A,(/)At )} 

+ /?3{Tr((/.A,0tAt) + Tr((/)tA,Mt)} . (8) 

All the parameters except a2 in the above are required to be real by imposing the discrete symmetry 

A^<->A„, ^^4>\ (9) 

simultaneously with the exchange of left-handed and right-handed fermions. Finally, a2 is chosen to be real from the 
requirement of spontaneous CP violation H^, |2|] . 

The ansatz for the VEV's of the scalar fields has been discussed extensively in the literature. After accounting for 
phases that can be eliminated by global symmetries and field redefinitions | |2C| ], only two independent phases remain. 
We choose them for convenience as follows in the translation invariant VEV's: 

where all the other parameters are taken to be real. Phenomenologically the hierarchy Vz, ki, k2, <^ is required. 
This separates the electroweak scale from the L-R symmetry breaking scale. It has been argued that this is 
possible to achieve for natural values of the above parameter set while also obtaining (1) spontaneous CP violation, 
(2) mixing of Wn, Zp^ with their SU (2) ^ counterparts which is unobservable at accessible energies, and (3) suppression 
of flavor changing neutral currents. 

Fermion masses are obtained from Yukawa couplings of quarks and leptons with the Higgs bosons. For one generation 
of quarks q and leptons ip, the couplings are given by 

Cy = h''qL(t)qR + h'' qL 4> qp. 

+ h'' IpLf/^TpR + h^ ^l4>^r 

+ f{^blC~'T2AL^bL+^lC-'T2AR^R)+h.C. (11) 

where C is the charge-conjugation matrix {ij}'^ = C'ljP"). Neutrino mass terms resulting from the above parameterization 
of the scalar VEV's are 

/:.-„ass = (/I'fcie*" + h'k2)v^, + {fvla''vpe'%^ + f^Wv^Vn + h.c] (12) 

The Majorana mass terms allowed for the neutrinos are a source of lepton number violation, while the CP violation 
needed for leptogenesis results from the phase a in the Dirac mass term. 



III. LEPTOGENESIS MECHANISM 



Sources of CP violation as well as existence of out-of-equilibrium conditions have been major challenges for realizing 
low energy baryogenesis. The presence of moving topological defects in unified theories is a novel source of out-of- 
equilibrium conditions. In ref. |lO[ it was shown that in the L-R model, at the first stage of gauge symmetry breaking, 
domain walls can form, which separate phases of broken 5C/(2)j^ and SU{2)l. The disappearance of the unstable 



5 



domains with unbroken SU{2)j^ provides a preferred direction for the motion of the domain walls. This can fulfill the 
out-of-thermal-equilibrium requirement for leptogenesis. 

Consider the interaction of neutrinos from the L-R wall, which is encroaching on the energetically disfavored phase. 
The left-handed neutrinos, Vl, are massive in this domain, whereas they are massless in the phase behind the wall. 
This can be seen from the Majorana mass term /i„Aii/f j/^, and the fact that (A^,) has a kink- like profile, being zero 
behind the wall and 0(wh) in front of it. 

To get leptogenesis, one needs an asymmetry in the reflection and transmission coefficients from the wall between 
Vl and its CP conjugate {vD- This can happen if a CP-violating condensate exists in the wall. This comes from 



the Dirac mass terms as discussed in section IV B. Then there will be a preference for transmission of, say, Vi,. 
The corresponding excess of antineutrinos {v'^) reflected in front of the wall will quickly equilibrate with due to 
helicity-flipping scatterings, whose amplitude is proportional to the large Majorana mass. However the transmitted 
excess of survives because it is not coupled to its CP conjugate in the region behind the wall, where (A^) = 0. 

A quantitative analysis of this effect can be made either in the framework of quantum mechanical reflection, valid 
for domain walls which are narrow compared to the particles' thermal de Broglie wavelengths, or using the classical 
force method |2^, which is gives the dominant contribution for walls with larger widths. We adopt the latter 
here. The thickness of the wall depends on the shape of the effective quartic potential and we shall here treat the case 
of thick walls. Further, we assume that the potential energy difference between the two kinds of vacua is small, for 
example suppressed by Planck scale effects. In this case the pressure difference across the phase boundary is expected 
to be small, leading to slowly moving walls. 

In refs. |2^, it is shown that the classical CP-violating force of the condensate on a fermion (in our case a 

neutrino) with momentum component perpendicular to the wall is 

F = ±sign{p,)^{mUx)a'{x))' . (13) 

The sign depends on whether the particle is ly^ or i/^, rn^{x) is the position-dependent mass, E the energy and a is 
the spatially varying CP-violating phase. One can then derive a diffusion equation for the chemical potential /x^ of 
the z^i, as seen in the wall rest frame: 

- D,^i'[ - v^fi'^ + e{x) Fm fi^ = S{x). (14) 

Here is the neutrino diffusion coefficient, is the velocity of the wall, taken to be moving in the +x direction, 
Fhf is the rate of helicity flipping interactions taking place in front of the wall (hence the step function 0{x)), and S 
is the source term, given by 

Six) = -^{y^Fix))', (15) 

where v is the neutrino velocity and the angular brackets indicate thermal averages. The net lepton number excess 
can then be calculated from the chemical potential resulting as the solution of eq. (|lj) . 

In order to use this formalism it is necessary to establish the presence of a position-dependent phase a. This is 
what we turn to in the following discussion of the nature of domain walls in the L-R model. 



IV. DOMAIN WALLS 



A. The Left-Right breaking phase transition 

The fundamental L-R symmetry of the model, eq. (^, implies that the gauge forces visible at low energies might 
have been the SU{2)ji rather than the SU{2)l with corresponding different hypercharge remnant of the U{1)b-l- In 
the early Universe when the symmetry breaking is first signaled, either A^, or A^ could acquire a VEV. In mutually 
uncorrelated horizon volumes, this choice is random. As such we expect a domain structure with either of these fields 
possessing a VEV in each domain. These may be referred to as L-like if they lead to observed phenomenology (with 
V-A currents), and R-like, if Aj, has remained zero. Such domains will be separated by domain walls, dubbed L-R 
walls in |l§. 
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The walls must disappear; otherwise they would contradict standard cosmology by dominating the energy density 
very soon after their formation. This must occur in such a way as to eliminate the i?-like regions. What biases the 
survival of the L-like regions cannot be predicted within the model. We will assume that there are small corrections 
suppressed by a grand unification scale mass, which favour the L-like regions slightly. A time asymmetry, due to the 
motion of the walls into the i?-like regions, arises as a result. The L-R walls necessarily convert the i?-like regions 
into L-like ones and disappear by mutual collisions. 

This can get implemented in two ways. One is explicit deviations from exact symmetry in the tree level Higgs 
lagrangian. An alternative is that the gauge couplings of the two SU {2)'s are not identical. In this case the thermal 
perturbative corrections to the Higgs field free energy will not be symmetric and the domain walls will be unstable. 

A possible reason for such small deviations from exact discrete symmetry could be that the model is actually 
descended from another unified model, and the small departures from exact symmetry are due to terms suppressed by 
the ratio of L-R breaking scale to the scale of higher unification. If the higher unification is in a conventional gauge 
group like SO{10), it may not constitute a good explanation since the breaking of such symmetry groups docs not 
generically result in a low energy model with close-to-exact L-R symmetry. It is however possible that the unification 
is of a different type, for instance supergravity or string unification, wherein mechanisms as yet not understood impose 
the kind of symmetry required, while permitting small energy differences in the free energies of the L-like versus i?-like 
phases. A study of disappearance of domain walls in the context of a Supersymmetric model has been made in 
and a study of the effectiveness of the mechancism in p6| . 

The breakdown of the L-R symmetry is described by the VEVs of two scalars A^,, A^. The form of the potential 
-(^ has been shown to have generic zero temperature vacua which are either i?-like or L-like. Let the difference in 
vacuum energy densities due to departure from exact L-R symmetry be SUl-r such that L-like vacuum is favoured. 
If this difference is purely in the scalar self couplings, it is determined directly by the GUT scale mechanism and will 
not be altered at finite temperature. On the other hand if the gauge couplings differ due to these suppressed GUT 
effects, the corresponding thermal corrections will thereby acquire differences, producing a corrected i5?7j_^ at finite 
temperature. The condition for the formation of the unstable domains can now be obtained as follows. If the phase 
transition is second order, its dynamics may be considered to have terminated after the Ginzburg temperature Tg is 
reached, which is given by ||2^ 

- A (16) 

where is the critical temperature and A the effective quartic coupling. The correlation length at this temperature 
is estimated to be — l/l-ATc). For the walls to form, the fluctuations that can convert the false vacuum to the 
true one must be suppressed before the Ginzburg temperature is reached. Thus the energy excess available within a 
correlation volume must be substantially less than the energy needed to overcome the barrier set by Tc — Tg, i.e., 

5ULaiG « T,-Tg (17) 

or 

5Ul_R « A^T^ ~ X^Vt^ (18) 

where we took the temperture-dependent VEV Vr to be ~ A^/^Tc. This bound is easily satisfied if the GUT scale is 
much higher than the L — R scale, as is expected. 



B. Wall profiles and CP violating condensate 

In order for nontrivial effects to be mediated by the walls, the fermion species of interest should get a space- 
dependent mass from the wall. Furthermore, the CP-violating phase a should also possess a nonvanishing gradient 
in the wall interior. We study the minimization of the total energy functional of the scalar sector with this in mind. 

The minimization conditions for the various VEVs introduced above are given in Appendix A, assuming trans- 
lational invariance. The presence of walls breaks this invariance, requiring derivative terms to be added in the 
minimization conditions. 
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We demonstrate that there are sizeable domains in the parameter space for which a position-dependent, CP-violating 
condensate results. In order to simplify the analysis we assume ki = k2 = k. The range of the parameter values for 
which such minima would be phenomenologically viable have been studied, e.g., in [ pT| . The analysis can be repeated 
for other cases along similar lines. Let the L-R wall be located in the y-z plane at a; = 0. Its equation of motion is 



da 



2 



fc-fc"-f l^^j k + {-nl-2filcosa + Afil)k 

+ {a2 + '^a3 + ^ai){vl+vl)k 

+ cos(q; — 6*) + /?3 cos 9 + 132 cos(2q; — 9)}vi,Vftk 

+ {Ai + A2(l + 4cos2q!) + A3 + 4A4Cosa}fc^ = . (19) 

The temperature correction to the mass-squared term (A/i^) is displayed explicitly. The remaining parameters are 
also mildly temperature-dependent but this is a small effect. The background fields w^, have solutions of the form 
At{1 ± tanh(x/A„)), with upper and lower signs being for L and R respectively. At is the temperature-dependent 
VEV which is possible for either of L or R fields. This value is of the order of the temperature T relevant to the epoch 
immediately after the L-R breaking phase transition. A„ is the wall width, of the order A;^^A~^/^, A here standing 
for the generic quartic coupling in the effective hamiltonian for the and fields. The nonderivative terms of this 
equation can be schematically written as 

m^k + A{vl+vl)k + Bv^vj,k + Lk^ = 0. (20) 

We are assuming > so that at the epoch in question, fc = in the absence of walls. This potential has two 
minima, 

k = or k = k(^o)^ -Y {m'^ + +vl) + Bvr^Vn) . (21) 

We want k ^ a,t the origin and fc = asymptotically. The latter is achieved if 

= AA^ > . (22) 

fe=0 



9fc2 

At the origin the nontrivial value fc(o) becomes 

fc(o)^"-^("^' + (2A + SKo)) ' (23) 
where W(o) — ^A-p is the common value of 1;^, at the origin. Thus 



k—k(i 



2Lfc^o) = -2(m2 + (2A + B)i;^o)) > 0. (24) 



Comparing eqs. ( P4 ) and (p4|), both conditions are satisfied provided the effective coefficient B becomes sufficiently 
negative. 

We can now proceed to determine a sufficient condition for a position-dependent nontrivial solution. We have 
already restricted ourselves to the case |fci| = fc2. We assume that the fates of the separate parts Im(fci) and Re(fci) 
are the same, i.e., if one of them is nontrivial, both would be so. So we focus on the condition for fc to be nontrivial. If 
the nontrivial solution is energetically favorable, the trivial solution should be unstable. Thus consider the linearized 
equation for the fluctuation 5k about the solution fc = 0. The desired time dependence of the solution is 

5k ~ 6^*=* X (spatial part) (25) 

with real parameter > for instability of the fluctuation. Then 

- 5k" + (m^ + A{vl + vl) + Bvr^Vn) 5k = -el 5k (26) 
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We compare this with the Schrodinger equation for a bound state wavefunction 

- ^"{x) + V{x)^p{x) = E^{x) (27) 

Our V{x) has the form of an attractive potential; it approaches a positive constant value Vq asymptotically, and 
Vo < near the origin due to eq. (|2^). In the Schrodinger equation above, for a bound state, £' < if V{x) 
asymptotically. In the present case, due to the positive constant value of V{x) asymptotically, bound states may exist 
even for E > 0. However our stability analysis demands E < 0, since we want efc to be real. If we ensure that the 
£" ~ solution has at least one node then there will be a lower energy solution with no nodes, the required unstable 
fluctuation. In the WKB approximation this condition amounts to 

^/^V{^)dx > - (28) 

where a and b are the zeros of V{x). Eqs. (p2|), ( ]2^ ) and ( ]2^ ) constitute one set of sufficient conditions on the 
parameter space for a CP violating condensate to occur within the width of the domain wall. They provide the range 
to be explored if a full numerical solution were to be attempted. 



V. EFFECTIVE HAMILTONIAN 



In this section we numerically study an effective Hamiltonian for the likelihood of generating a CP violating 
condensate. In a suggestive notation we choose three fields L, R and K representing the VEVs of A^,, A^. and the 
electroweak Higgs respectively. The energy per unit area of the wall configuration can be taken to be 

H = ldx{^\Lf + ^p,\L\'{\L\-Mf + ^{Rr + lpiR'{R-Mf+p,\L\'R' (29) 

+ + ^X{\K\^ + m^f + ailif p(|L|2 + R^) + p^\K\''{ReL)R 

+l32\K\Re{KL)R + fi3Re{K^L)R} 

where L and K are complex and R is real. M represents the Left-Right breaking mass scale and m the electroweak 
breaking mass scale, both including the appropriate finite temperature corrections. Thus rn? is positive. Likewise the 
other parameters are determined by the parameters of the original lagrangian. The equations we get after rescaling 
the fields by M are 



PiLl 1 - 



2y/ Lj + L 



+Li 



-R 



Pi 



^2 



Ll 



^2 



2p;R^ 



/3i [Ki + Ki) + P^K^d Kl 



2yZfTX|y 



2a (Kf + Kl) 



dx^ 



PiLl 1 



Ll. 



-R 



Pi 



2^ 



Ll- 



Ll 



Ll, 



2p^R^ + 2a {Kf + Kf) 



f32K2^Kl + Kl + 2/33i^ii^2 



0.8 - 



§0.6 

> 

-a 

I 0.4 



0.2 



Q 
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X = 0.9 
a = 0.l 

Pl=-0.9 
P2=-0.5 
P3=0.5 
Pi=0.1 
P3=0.5 
m/M = 0.01 



FIG. 1: L, R and K values for given choice of ps, ai and (3i, in M units, x values in M ^ units. Values of parameters not 
shown are unity. 



= R 



pii? ( i? - - 1 + ^ + 2p3 (i? + Ll) + 2a {Kf + Kl) 



l3iLi {Kl + Kl) + 132 J Kl + Kl [K^Li - K2L2) 



dx^ 



X[Kl + Kl + j' ) + 2a {L\ + LI + R^)+2 (/3i + /^g) RL^ 



X[Kf + Ki+[^y]+2a{Li + Ll + R^)+2{Pi-Ps)LiR 



+K,MEl^^^ - P2RL2J^^[^l^2p,RL2Ki 




In addition to the above we need the expression 

-7—7 = 2ai(|L|2 + i?2) + 2/3ii?ii + 2/32i?(Licosa-L2sina) (30) 
a\K Y 

+2P3R (Li cos 2a - L2 sin 2a) + X{3K^ + rri^) 
for studying stabihty issues. This shows that to ensure K=0 asymptotically (no E-W breaking at L-R scale) we need 

2«, + a(-) >0. 

To obtain X 7^ in the core of the wall, we again use (|3^) with the values of R and Li in the core estimated to be 
0.5. This suggests the requirement 

(Til 1 
^) --(2ai+/3i+/32 + /33)>0 

This has to be revised in view of the actual values of R and Li due to backreaction of the K fields, but serves as a 
good indicator to the required range of values. 

The second equation clearly suggests taking /3's negative. In particular (3i can be 0(1) and negative, which will 
ensure the required instability of the K = Q vacuum inside the wall core. Two examples of numerical solutions are 
shown in figures |l| and Parameters other than those displayed are taken to be unity. The shape of the K profiles 



11 



is of the form of the sech function, as expected for lowest hnear perturbation in a tanh background. The numerical 
study indicates the field profiles are sensitive to the parameters governing the Yukawa couplings, but they have no 
appreciable variation with respect to the ratio of the mass scales (■^) . 



VI. IMPLICATIONS FOR COSMOLOGY 



We are now in a position to use the formalism of eqs (|T^)-(|T^) to estimate the lepton asymmetry generated. The 
asymmetry in local number density is given by 



6 

where /i^ satisfies the diffusion equation ( [T^ ) . The general form of the solution to this equation is 



p+(x-y) _ pP- i'-c-y) ' 



Siy), x>0 



where 



B = 



0, A_ 
1 



Em. 

D 



D{\--K) 



A = 



1 



-1/2 



(31) 
(32) 

(33) 
(34) 



The integration constants A^ and B^ are chosen so that pix) and its derivative are continuous at a; = 0, and /j, is 
finite as a; — > ±cx3. In particular, we are interested in the limiting value /xq = lima;^_oo /it(x) = since this is 
relevant deep within the L-like phase and represents the uniform lepton asymmetry filling the universe long after the 
wall has passed by. It is given by 



1 



Mo 



dy 1 



P+ 
— I 

P- 



S{y) 



1 



P-D 



''Sl^y). 



(35) 



We note that in the limit Vw 0, the above expression is finite for a generic source S{y). But since our source 
vanishes when Vw — 0, we get no lepton asymmetry in that limit, which is in accord with Sakharov's out-of-equilibrium 
requirement. We can also verify that no lepton asymmetry arises when lepton violating interactions are turned off. 
For us, this means setting Fhf = 0, in which case we obtain /io — S{y)dy. The integral vanishes if the source 

itself does not violate global lepton number conservation, so this check also succeeds. The third necessary ingredient, 
CP violation, is contained within the source Sly), since this depends on the neutrino masses having complex phases 
which vary within the domain wall. 

Now we proceed to estimate the chemical potential po which quantifies the generated lepton asymmetry. This 
requires the thermal averages |2q| 



3xi/ 



£2 



X ^ 3x u 

a — bxi, 
2^2 ' 



- 1; 



-2 
a ^0.24: 



— 

6 = 0.65 



(36) 
(37) 



The first approximation (^) is good for relativistic neutrinos with x^, <, 0.1, and the second one ( pTj ) is an approx- 
imation to the function given in [^8| which is adequate for our estimate. By taking {v^) ^ 1 we can simplify the 
expression for since the source S{y) becomes a total derivative. Integrating by parts, 



Mo 



Vw P+ 

r2 p_ 



-P+V 



rp2 



{m^a'Y dy 



(38) 



12 



Since the wall is much thinner than the diffusion scales D jv^ and l/p+ , it is a good approximation to neglect these in 
the integral (setting e^^^l^ and e~P+y to 1). We will use the ansatz m^^y) = M^h'^iy), h{y) = i(l + tanh(?//A^,)), for 
the real part of the neutrino mass, while for the phase, in accordance with the profiles found in the previous section, we 
take a{y) = Im(L(y))/Re(L(t/)), with Im(L(y)) = q:oA^,/i'(j/). Here is the large value of the neutrino Majorana 
mass neutrino, acquired by the left (right) -handed neutrino in the i?-like (L-like) phase. We have performed the 
integral numerically to obtain the analytic result: 

^n,^Omv^^^ (39) 

This is the raw value of the lepton number generated by this mechanism. One would like to express this as a ratio r/^, 
of lepton number to entropy as is standard to do with baryon number. Using the expression for the entropy density 
S = 2tt'^ g„T^ / 45 of g„ relativistic degrees of freedom, we get 

=0.01z;„^-^ (40) 

Let us consider whether this result can naturally be of the same order as the observed baryon asymmetry. Let Aj^ 
[as introduced above, eq. (^9|)] denote the temperature-dependent VEV acquired by the A^ in the phase of interest. 
Experience with the electroweak theory shows that Aj^/T is determined by the ratio of gauge and Higgs couplings, 
and is typically smaller than unity. If / is the Yukawa coupling determining the Majorana mass, then M^, — /At- 
Moreover the inverse wall width A^^ is V^^cff At, where Aoff is the effective quartic self-couplirig of the A fields. This 
is assumed to be small, since we have taken the wall to be thick. Therefore we can reexpress ( pOj ) as 

^raw 




ryr = 0.01z;„,- -f fWKs- (41) 



With Ki 10^, this raw lepton asymmetry is close to rj^ = 10 the desirable value for final baryon asymmetry, 
provided that 

(^)'/^Va;J^-10-. (42) 

Even if At/T ~ 0(1) and \/Aoff ^ 1, this can be achieved with a reasonable value for the Majorana Yukawa coupling 
/a « 10~^ of the heaviest (third generation) sterile neutrino. Ignoring further evolution of the lepton asymmetry for 
the moment, one could turn this around to derive a lower bound on /, assuming that the present mechanism was 
responsible for baryogencsis. If all the Majorana neutrinos are lighter than « IO^^vr, then it produces too small a 
lepton asymmetry to be significant. 

After the domain walls have disappeared, the lepton asymmetry undergoes further processing by several interactions. 
Firstly the electroweak sphalerons will redistribute this asymmetry partially into baryonic form. This is the mechanism 
by which we get baryon asymmetry from the wall generated lepton asymmetry. The standard chemical equilibrium 
calculation ||2^ gives 

28 28 
Ann = —AuB-L = -TT^n^ (43) 
79 51 

assuming the minimal Higgs and flavour content of the Standard Model. 

However the presence of heavy majorana neutrinos gives rise to processes that can deplete the lepton asymmetry 
generated. Such processes were considered in a model independent way in jl^, |l5|, referred to in section I. The 
present model differs from classic GUT scenarios in that the temperature Tb-l = ^lr when the lepton asymmetry 
is created can be less than or comparable to the heavy neutrino mass Mjv . The two processes of importance are the 
decay of the heavy neutrino with rate Fq and heavy neutrino-mediated scattering processes with rate Fg. The latter 
class of processes in the context of the present model is shown in fig. ^. The rates are given roughly by 

r " r ^ _! (AA) 

° 167r(4r2-HM2)i/2 ^ 137r3 (9r2 M2 ) ^ ^ 
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These expressions correctly interpolate between the high and low temperature limits which can be inferred from eqs. 
(3.1,3.8), (A15-A16) of ref. |jl^, using the Boltzmann approximation Ki{x) ^ jx in the thermal average of the 
scattering cross section. (The factor 13 in (^) is really 96C(3)/9.) 

Let us first consider the case when the decays do not deplete the generated lepton asymmetry at all. This happens 
if the lightest of the heavy majorana neutrinos has Mjvi > T^^h, so that the decays do not occur because of Boltzmann 
suppression. This limit tends to make the initial lepton asymmetry 77^, large, possibly 0(1) from (^0|). However the 

lepton-violating scattering processes will dilute this by the factor 10^''^ = cxp IT^r ^ sdt^ where tj^Yi is the time 
of the LR-breaking phase transition and is the present. At the same time, sphalerons will keep the baryon and 
lepton asymmetries in the same proportion p9[ until the electroweak phase transition, at which time the sphalerons 
go out of equilibrium. The corresponding depletion factor for baryons, rewritten in terms of an integral with respect 
to temperature, is 

C/'T'min p (iT\ 
- -j^jrj ; T,„i„ = min(r^H, r.ph) (45) 

where T^ph ~ 10^^ GeV is the maximum temperature below which sphalerons are in equilibrium. Evaluating the 
integral gives the baryon depletion exponent 

'^-^ 13.MnlO./a: ^^i^ (4') 



where is the average number of relativistic degrees of freedom, and we are assuming that Af„i > Tmin 

can be solved for the Yukawa coupling h which gives the Dirac mass term for the neutrino: /i^ <, 3200 dg i ^ j^j^ 

where we have taken g^, — 110 for definiteness. Since dg should be no greater than about 10 to avoid too much 
dilution of the baryon asymmetry, this can be further transformed into an upper limit on the light neutrino masses 
using the seesaw relation ~ (hv)'^ /Mn where v is the Higgs boson VEV, v = 174 GeV: 




™' ^ ^T~M^i ^^^^ 

If the heaviest neutrino mass is 1 eV, for example, the temperature of the LR phase transition (if it is smaller than 
Tsph), being also the temperature at which most of the _B — i is generated, is predicted to be 

T._. .r.„ <10.»G.Vx(£)'x(^£) (48, 




FIG. 3: The mediated processes violating L 

The previous discussion of dilution by lepton-violating scattering assumed the heavy neutrinos N had masses 
Mjv > so that the decay processes could be neglected. If we are in the opposite regime, < T^r, the decays 
and inverse decays of N will dominate over scattering for the epoch of temperatures T > Mjy. For lower temperatures. 
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the decay rate is exponentially suppressed by the Boltzmann factor e ^^«/^. We can roughly estimate the dilution 
due to decays as 




in the limit that Mjv ^ Timn- Again requiring that < 10 gives the bound on the heaviest neutrino mass 



■nil, < 0.3 eV x 




) 



(50) 



It is interesting that this value is compatible with, and not very far from the value implied by atmospheric neutrino 
observations. 



We have shown that a hitherto unexplored mechanism exists in the Left-Right symmetric model for generation of the 
observed baryon asymmetry of the Universe. The idea is reminiscent of electroweak baryogenesis, but here the motion 
of domain walls with CP-violating reflections of neutrinos during the LR-breaking phase transition creates a large 
lepton asymmetry, which is subsequently reprocessed by sphalcrons into the baryon asymmetry. Unlike electroweak 
baryogenesis, there is no suppression by a^, since the sphalerons are in equilibrium and they have sufhcient time 
to equilibrate the baryon and lepton numbers. Rather, the answer is determined to a large extent by (see (|4l|)) 
because asymmetry production is determined by the helicity flipping interactions. There are no natural smallness 
requirements on this parameter, although through see-saw formula it is constrained by the observed light neutrino 
mass. Furthermore there are no strong constraints on the CP violating phases since they appear in the interactions 
of the heavy right handed neutrino. 

It is possible to generate the observed baryon asymmetry for a range of parameters of the model. We have studied a 
few limiting cases to demonstrate the intrinsic potential of this scenario for producing the observed baryon asymmetry. 
One extreme possibility is that we could start with the raw value of the lepton asymmetry ( p] ) being of order 10" 
by virtue of a small Majorana Yukawa coupling, eq. (^2|), while the heaviest left-handed neutrino mass satisfles the 
bounds (^, |5^) (evaluated at ds ~ 1) which guarantee that there is no subsequent dilution of the asymmetry by 
Icpton-violating interactions. The other limiting case is to initially create an asymmetry of 0(1), by taking large 
Majorana Yukawa couplings; the asymmetry is subsequently diluted to the required level by saturating the bounds 
( ^ , |50| ), which make reference to the heaviest left-handed neutrino mass. 

An interesting application of this mechanism is the possibility to generate a large lepton number as suggested in 
pO| and considered in the context of new observations in cosmology e.g. in [|T]-|3^, notably the MAP and PLANCK 
experiments to measure the cosmic microwave background (CMB) fluctuations. In the simplest model with just one 
lepton generation, we cannot create a large lepton asymmetry without also making the baryon asymmetry too large. 
But consider a model with a certain combination of lepton numbers conserved, such as Lf.-\-L^. This would be the case 
if the Majorana mass matrix had the form (^^q). Then the leptogenesis mechanism would create equal and opposite 
amounts of and L^. Since sphalerons separately conserve — and — L^, the combination ^B — — 
would remain conserved at all times, so that the resulting baryon asymmetry would be zero even if \Le\ and \L^\ 
separately were large. By adding a very small breaking of the + symmetry, one could generate the observed 
baryon asymmetry simlutaneously with large lepton asymmetries In addition to its imprint on the CMB, such 
an effect could have other observable consequences as observations relevant to nucleosynthesis are improved - [Eoj . 
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Appendix A: minimization conditions for wall profiles 



These conditions for finding the wall profiles were used in section IV B 
5V 



Ski 



5V_ 

5a 



+ 2piw^ + pst'twi + 4Q!2fcifc2Wt cos(q!) 
+ [3ikik2VRCOs{a — 6) + f32k'lvaCOB{2a - 

= ]j^{oiiklvji + aiklvi,+ a^iklv,^ 

+ 402^1^2^^ cos(a) + /3ifcifc2fi cos(a — 6) 
+ (32klvLCOs{2a- e) + (33klvi^cos{9)} 

^ = -^v lV „{- I3ikik2 sm{a - e) - l32kl sm{2a - 9) + Pskj sin 0} 
00 2 

5V 

— = Xi{kfk2 + k^) + 2X3klk2 
ok2 

^i\k2 + ^aik2vl + ^a3k2vl 
iaifc2W^ + ^ai.k2vl, + k^X^ cos a 

ikik^Xii cos a — 2kiii\ cos a + a2kiv'^ cos a + a2kiv\ cos a 



Svt 



SV 



+ 4fcjfc2^a2 cos 2a + -Pikiv^Vji cos(a — 6) + /33A:2^'i,^'r cos ( 



— = k'fXi + kikjXi + 2kiklX3 - kifif 



—aikiv^ + —aikiv^ + 3kfk2X4 cos(q!) 

A4 cos(q!) — 2fc2/i2 cos(a) + 02^2'^'^ cos(q!) 
(^2^2^^ cos(a) + 4fciA:2 A2 cos(2q;) + ^/3ifc2Wi,UH cos(q; — 9) 



P2kiVLVii cos(2q; — 9) 

— A4/cifc2(fc^ + /cj) sina + 2kik2ii\ sina 

o-2kik2{v\ + v\) sin(a) — 4fc^fc2A2 sin 2a — —(3ikik2VLVn sin{a — 9) 
/32^i Wt^H sin(2Q! — 9) 
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